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We define a cellular automaton where a resting cell excites if number of its excited neighbours 
belong to some specified interval and boundaries of the interval change depending on ratio of excited 
and refractory neighbours in the cell's neighbourhood. We calculate excitability of a cell as a number 
of possible neighbourhood configurations that excite the resting cell. We call cells with maximal 
values of excitability conductive. In exhaustive search of functions of excitation interval updates 
we select functions which lead to formation of connected configurations of conductive cells. The 
functions discovered are used to design conductive, wire-like, pathways in initially non-conductive 
arrays of cells. We demonstrate that by positioning seeds of growing conductive pathways it is 
possible to implement a wide range of routing operations, including reflection of wires, stopping 
wires, formation of conductive bridges and generation of new wires in the result of collision. The 
findings presented may be applied in designing conductive circuits in excitable non-linear media, 
reaction-diffusion chemical systems, neural tissue and assembles of conductive polymers. 

I. INTRODUCTION 

Excitable cellular automata are well endowed tools for studying complex phenomena of spatio-temporal physical, 
chemical and biological systems [U |2] , prototyping of chemical media [3l |4] , react ion- diffusion computers [5] , studying 
calcium wave dynamics [6 , and chemical turbulence [7 . 

In a classical Greenberg-Hasting [8 automaton model of excitation a cell takes three states — reseting, excited and 
refractory. A resting cell becomes excited if number of excited neighbours exceeds a certain threshold, an excited 
cell becomes refractory, and a refractory cell returns to its original resting state. In [9 we introduced a bit more 
exotic cellular automaton, where a resting cell is excited if a number of its excited neighbours belongs to some fixed 
interval [6>i,^2]- The interval [^1,6^2] is called an excitation interval. For a two-dimensional cellular automaton with 
eight-cell neighbourhood boundaries of the excitation interval are between 1 and 8: 1 < ^1 < ^2 ^ 8- By tuning Oi 
and O2 we control automata dynamics and evoke target and spiral waves, stationary excitation patterns, and mobile 
localisaitons [5 . 

How does excitation dynamics change if we allow boundaries of the excitation interval to change during the au- 
tomaton development? We partially answered the question in [10 by making the interval [0\{x) , OHx)] of every cell 
X dynamically updatable at every step t depending on state of the cell x and numbers of excited and refractory 
neighbours in the cell x's neighbourhood. We found that excitable cellular automata with dynamical excitation in- 
tervals exhibit a wide range of space-time dynamics based on an interplay between propagating excitation patterns 
and excitability of cells modified by the excitation patterns. Such interactions lead to formation of standing domains 
of excitation, stationary waves and localised excitations. We analysed morphological and generative diversities of the 
functions studied and characterised the functions with highest values of the diversities. 

Excitable cellular automata with dynamical intervals of excitation can be considered as discrete phenomenological 
models, or rather conceptual analogs, of memristive media and excitable chemical medium computers. 

A. Memristive medium 

The memristor — a passive resistor with memory — is a device whose resistance changes depending on the polarity 
and magnitude of a voltage applied to the device's terminals and the duration of this voltage's application. Its existence 
was theoretically postulated by Leon Chua in 1971 based on symmetry in integral variations of Ohms laws |1THT3] . 
The first experimental prototypes of memristors are reported in [ 14HT6] . An importance, and the great pragmatic 
value, of memristors is that one can design logically universal, or functionally complete, circuits composed entirely of 
the memristors. Potential unique applications of memristors are in spintronic devices, ultra-dense information storage, 
neuro-morphic circuits [17 , and programmable electronics [18 , designing binary arithmetical circuits with polymer 
organic memristors [19]. 

Despite phenomenal number of results in memristors produced literally every week there are insufficient findings on 
phenomenology of spatially extended non-linear media with hundreds of thousands of locally connected memristive 
elements. Three cellular automaton models of a memristive medium were suggested so far: 

• Itoh-Chua memristor cellular automata, where cellular automaton lattice is actually designed of memristors [20] , 
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• Adamatzky- Chua model of memrisitive cellular automata based on structurally-dynamic cellular automata [21], 

• semi-memristive automata p2] . 

Itoh-Chua and Adamatzky- Chua models imitate memristive properties of links, connections between cells of au- 
tomata arrays but not the cells themselves. The semi-memristive automata bring memristivity into cells [22 : links 
between cells are always 'conductive' but cells themselves can take non-conductive, or refractory, states. The semi- 
memristive automata are excitable cellular automata with retained refractoriness. 

In present paper, an excitability of a cell is calculated as a number of possible neighbourhood configurations that 
excite the resting cell. Cells with maximal values of excitability are called conductive. We represent conductivity 
of a cell X via boundaries of its excitation interval [0i{x)^02{x)]. We excite the cellular automaton, wait till the 
perturbation spreads and boundaries of excitation intervals of cells updated and then select cells with highest values 
of excitability. The configurations of conductive cells form conductive pathways by analogy with the formation of 
conductive pathways in disordered networks of organic memristors [23 . In automata studied polarity of a direct 
current applied to a cell x is imitated by excited and refractory neighbours of x, and current intensity is represented 
by a ratio of excited and refractory neighbours. 



B. Ensembles of Belousov-Zhabotinsky vesicles 

Excitable reaction-diffusion computers, especially those based on Belousov-Zhabotinsky reaction, employ principles 
of a collision-based computing [5 . Wave-fragments collide in a 'free' space and change their velocity vectors in the 
result the collision. When input and output waves are interpreted as logical variables., the site of the waves' collision 
can be seen as a logical gates. Wave-fragments, similar to dissipative solitons [24 , are inherently unstable: they 
either collapse or explode. A way to overcome the problem of wave-fragments' instability was suggested in [25H27] : a 
subdivision of the computing substrate into interconnected compartments, so called BZ-vesicles, and allowing waves 
to collide only inside the compartments. Each BZ-vesicle has a membrane that is impassable for excitation. A pore, 
or a channel, between two vesicles is formed when two vesicles come into direct contact. The pore is small such that 
when a wave passes through the pore there is insufficient time for the wave to expand or collapse before interacting 
with other waves entering through adjacent pores, or sites of contact. 

It has been observed in chemical laboratory experiments with BZ-vesicles [28] that waves of oxidation induced by 
external stimulation, e.g. with a silver wire, propagating in the initially resting BZ medium may cause changes in 
the excitability of BZ-vesicles (not just refractoriness but excitability in a long run). Excitability of BZ-vesicles can 
increase after first round of the oxidation wave propagation. A cellular automaton model designed in present paper 
gives a phenomenological snapshot of an ensemble of regularly arranged BZ-vesicles (imitated by cells), which change 
their long-term excitability after being subjected to propagating waves of excitation. 

The paper is structured as follows. We define an excitable cellular automata with dynamically updated boundaries of 



excitation interval in Sect.jllj Configurations of conductivity generated by the automata are analysed in Sect.|III[ The 
functions which produce fully conductive configurations are selected in exhaustive search. Section [IV| demonstrates 
how to design conductive wire-like pathways by positions seeds of excitation. Potential further developments are 
outlined in Sect. IVl 



II. EXCITATION CONTROLLED EXCITATION INTERVALS 

Let L be a two-dimensional orthogonal array of finite state machines, or cells, x'^ and x*~^^ states of a cell x at time 
steps t and t + and cr^(x) a sum of excited neighbours in cell x's eight-cell neighbourhood u{x) = {y : \x — y\L^ = 1}. 
Cell X updates its state by the following rule, = f{u{x)^ where cell-state update function / is represented as 



, ifx^ = o and cr^(x)+ G [e{{x), 6>^(x)] 

, ff = + (1) 
otherwise 



A resting cell is excited if a number of its neighbours belongs to excitation interval [^i(x), ^2(^)]? where 1 < 
6{{x)^62{x) < 8. The boundaries 6l{x) and OKx) are dynamically updated depending on cell x's state and num- 
bers of x's excited cr^(a;) and refractory <jL(x) neighbours. A natural way to update boundaries is by increasing or 
decreasing their values as follows: 

0{+\x) = mi^) + Ai(/.K(x) - aUx))) 
0'+\x) = mix) + A2^K(x) - al{x))) 
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where 




T2, if x = - 
T4, ifx = - 
0, if X = 



(3) 



and (j){a — 6) = 1 if a > 6, if a = 6 and -1 if a < 6, and ^(a) = 1 if a < 1 and 8 if a > 8. Boundaries of excitation 
interval [0{{x)^ ^ii^)] updated independently of each other. Rules of excitation intervals update are determined by 
values of Ti, • • • T4. We therefore address the functions as tuples £^(Ti, T2, T3, T4) which range from £^(—1,-1,— 1,-1) 
to 1,1,1). 

Excitability f (6>i(x), 6>2(x)) of a cell x with excitation interval [0i{x)^62{x)] is measured as a number of all possible 
local configurations which have sum of excited cells lying in the excitation interval [0i{x)^02{x)]: 



m, ^2) = |> e {o, +,-}«}: /M = +1 



(4) 



Two highest excitability values are reached by cells with excitation intervals [1,7] and [1,8], f(l,7) = 6304 and 
f (1,8) = 6305. We assume a cell x is conductive if £{6i{x)^62{x)) > 6300. Connectivity of the configurations of cells 
with maximal excitability is considered to be analog of conductivity of the whole cellular array. The connectivity is 
estimated using standard bucket fill algorithm. 

Definition 1 We call conductivity configuration diameter D fully conductive if there is a path between two sites lying 
at distance D from each other, and over 90% of sites are reached from any given site. 

We study formation of conductive pathways in initially non-conductive medium. Therefore, in experiments we 
considered initial conditions Oi{x) = 2 and ^2(^) — ^ for any x. The excitation interval [2,8] gives a cell excitability 
5281, which is below a threshold adopted as indicator of conductivity. 

We excite a disc radius 200 cells with a random configuration of excited states. Let p be a probability for a cell x 
to be assigned excited state at t = and one of its neighbours, chosen at random, also assigned an excited state. Two 
neighbouring excited cells is a minimal size of perturbation due to O^^x) — 2. We considered p = 10~^ and p = 0.1. 
Initially cells inside a disc radius 200 are assigned excited states with probability p. 

In each trial we allowed the cellular automaton to develop for 440 iterations and then analysed configurations of 
conductivity. Size of cellular arrays was chosen large enough for a front of perturbation to never reach the array's 
boundaries in 440 steps, so there are no infiuences of boundary conditions. 



III. CONDUCTIVE CONFIGURATIONS 



Sites of initial random excitations generate waves, localizations and other travelling patterns of excitation. The 
patterns either stay localised or merge and propagate outwards the initially perturbed region. The patterns of 
excitation update boundaries of excitation intervals of cells they excited. An example of cellular automaton, function 
£'(— 1, 1,0, —1), development is shown in Fig. [l] Waves of excitation propagate from the perturbation sites (Figs, [l^ 
and[2|, leaving somewhat fibre-like trails and extended domains of lower 6>i (Fig.[l]3) and upper (Fig.jl]^) boundaries 
of excitation interval. These are reflected in solid domains of conductivity partially linked with wire-like conductive 
paths (Fig.[l]i). Excitation waves originated from different sites of perturbation merge outside the stimulation disc and 
propagate further as almost connected packet of target waves (Fig.[l^). These waves leave triangular solid domains 
of conductivity behind (Fig. [l]i). 



A. Classes of connectivity 

We characterise local connectivity of conductivity configurations using z^max, a number of conductive neighbours 
of a conductive cell occurred most frequently in the configuration of conductivity, and z^min, occurred less frequently. 
For initial probability of excitation 10~^ we have twelve classes of local connectivity, the functions are grouped by 
values M = (z^max, ^min) of their conductivity configurations. Examples of configurations of conductivity generated 
by functions from these classes are shown in Figs. |3] and [4] 



• AT = (0,0): ^(0,-1,0,-1), ^(0,-1,0,0), ^(0,-1,0,1), ^(0,0,0,-1), ^(0,0,0,0), ^(0,0,0,1), E(0, 1,0,-1), 
E(0, 1,0,0), ^(0,1,0,1). 
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(d) Conductivity 

FIG. 1: (Color online) Snapshots of (a) excitation pattern, (b) configuration of lower boundary Oi of excitation interval, 
(c) configuration of upper boundary O2 of excitation interval, (d) configuration of conductivity of cellular automaton governed 
by function ^( — 1, 1,0, —1). Array of 880 x 880 cells evolved for 440 steps. Initially cells inside a disc radius 200 are assigned 
excited states with probability 10~^. Cell states in (a) are represented by colours and grey levels as follows: excited state + is 
red (c. 76 grey), resting state is white and refractory state — is blue (c. 28 grey). Colour values of excitation interval boundaries 
Oi (b) and O2 (c) are following: 1 is white, 2 is green or 150 grey, 3 is yellow or 226 grey, 4 is blue or 28 grey, 5 is magenta or 
104 grey, 6 is cyan or 178 grey, 7 is red or 76 grey, and 8 is black. Conductive cells in (d) are black and non-conductive are 
white. 

• A/- =(0,4): £;(1,0,0,-I),i3(l, 1,0,-1) (Fig.|§i) 

• A/" = (2, 5): E{<d, -1, 1, 1) and £^(1, 1, 0, 0) (Fig. ^) 

• A/- =(2, 6): i?(l, 1,0,1) (Fig.|§) 

. = (2, 7): E(0, 0, 1, 1), E(0, 1, 1, 1), E{1, -1, 0, -1), E{1, -1, 1, -1), E{1, -1, 1, 0), E{1, -1, 1, 1), E{1, 0, 1, -1), 
E(l, 0,1,0), E(l, 0,1,1), E(l, 1,1,-1), £(1,1,1,0), £(1,1,1,1) (Fig.|3}l) 

. AA = (2, 8): £(0, -1, 1, -1), £(0, -1, 1, 0), £(0, 0, 1, -1), £(0, 0, 1, 0), £(0, 1, 1, -1), £(0, 1, 1, 0) (Fig. ^) 

. AA=(5,0): £(1,-1,0,1), £(1,0,0,1) (Fig.jf) 

. AA=(7,0): £(1,-1,0,0), £(1,0,0,0) (Fig.|§i) 
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FIG. 2: (Color online) Enlarged part of excitation pattern from Fig. [l 



. Ar= (8,0): ^(-1,-1,-1,-1),^(-1,-1,1,1),E(-1,0,-1,-1),E(-1,0,0,-1),^(-1,0,0,0),E(-1, 1,0,-1), 
i?(-l, 1,0,0), i?(-l, 1,0,1), i?(0, -1,-1,-1), i?(0, -1,-1,0), i?(0, -1,-1,1), i?(0, 0,-1,0), i?(0, 0,-1,1), 
£;(0, 1,-1,1), £(1,-1,-1,-1), i?(l, -1,-1,0), i3(l, -1,-1,1), i?(l,0,-l,0), £;(1,1,-1,1) (Fig.^) 

. AA= (8,1): £(-1,-1, -1,0), £(-1,-1, -1,1), £(-1,0, -1,0), £(-1,0, -1,1), £(-1,1, -1,0), £(-1,1, -1,1), 
£(1,0,-1,1) (Fig. 15) 

• Af = (8,3): £(-1,-1,0,0), £(-1,-1,0,1), £(-1,-1,1,0), £(-1,0,1,1), £(-1,1,-1,-1), £(-1,1,1,1), 
£(0,0,-l,-l),£(0,l,-l,-l), £(0,1,-1,0), £(1,0,-1,-1), £(1,1,-1,-1), £(1,1,-1,0) (Fig.jiji) 

. = (8,6): £(-1,-1,0,-1), £(-1,-1,1,-1), £(-1,0,0,1), £(-1,0,1,-1), £(-1,0,1,0), £(-1,1,1,-1), 
£(-1,1,1,0) (Fig. 15) 

Class Af = (8, 0) is the largest class, it has 19 functions. Configurations of conductivity generated by the functions 
from Af = (8, 0) (Fig. ^p) are characterised by large solid conductive domains formed either by initial source of 
excitation, see e.g. rectangular embedded shapes in Fig. ^p, or by merging fronts of excitation (solid polygonal 
domains at the periphery). 

The next largest classes are Af = (2, 7) and Af = (8, 3) (Fig. |4ji)- Each of them includes 12 functions. Conductivity 
configurations generated by functions of class Af = (2, 7) consist mainly of isolated line segments of cells in conductive 
states and pairs or singletons. Configurations generated by functions of Af = (8, 3) are comprised of solid domains of 
conductive states and sparsely scattered singletons (Fig. |3ji) . 

Classes AT = (2, 5) (Fig.|3]D), A/" = (2,6) (Fig.|3];) andAr= (2,7) (Fig.jsji) pose particular interest because majority 
of conductive cells have two conductive neighbours each and therefore chances of conductive 'wires' to be formed 
could be high. 
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(d)^(l, 1, 1, 1), N = (2, 7) (e)^(O, -1, 1, -1), Af = (2, 8) (f)^(l, -1, 0, 1), Af = (5, 0) 



FIG. 3: Examples of conductivity configurations. 

Configurations generated by functions from J\f = (2, 5) (Fig. [sJd) predominantly consist of singletons and pairs of 
cells in conductive states; the pairs are aligned in staircase like fashion but there is no immediate connection between 
the pairs (Fig. [S^i). There are zig-zag like 'wires' of conductive states (Fig. [s})) in configurations (Fig. [sj^) generated 
by the only function ^(1, 1, 0, 1) in class JV = (2, 6). However they are present as isolated fragments. Configurations 
(Fig. [3^) generated by function £^(0, 1, 1,0) of class JV = (2,8) consist of clusters of singletons — single conductive 
cells surrounded by non-conductive neighbours — and pairs of cells in conductive states. The clusters are connected 
to each other with diagonals of singletons (Fig.lsfc). Neither of the functions with Vrnax = 2 generate fully conductive 
configurations but few functions from classes A/ = (7, 0) and A/" = (8, 0) do. We study these functions below. 

B. Functions generating conductive configurations 

Most functions generate only partially conductive configurations of conductivity. Examples of connected clusters 
are shown in Fig. |6] 

Finding 1 For low probability of initial excitation, p = 10"^^ functions £^(1,0,0,0) and £(1, —1,0,0) generate fully 
conductive configurations. For high probability of initial excitations, p = 0.1^ fully conductive configurations are 
generated by functions £(—1, 0, 0, 0) and £( — 1, 1, 0, 0). 

For low probability p = 10~^ of initial excitation only functions £(1, —1,0,0) (Figs. [7^) and £(1,0,0,0) generate 
fully conductive configurations. For initial configurations with high probability of excited cells {p = 0.1) only two 
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(a)£;(l, -1,0,0), AT = (7,0) 



(b)£;(l, -1,-1,1), AT = (8,0) 






(c)^(l,0,-l,l), A/'= (8,1) 



(d)£;(-l,-l,0,0), A/'= (8,3) 



(e)£;(-l,-l,0,-l), A/'= (8,6) 



FIG. 4: Examples of conductivity configurations. 

functions generate conductive configurations: £^(—1, 0, 0, 0) (Fig.jTj^) and £^( — 1, 1, 0, 0) (Fig. [7]i). In all four functions 
selected values T3 and T4 are nil. The values correspond to update of excitation intervals of refractory cells. This is 
in line with commonly accepted assumption that in an excitable medium elements in refractory state are insensitive 
to states of their neighbours. Mechanics of excitation interval updates of excited cells is illustrated in Fig. [Sj 

Function £^(1,0,0,0) increases lower boundary Oi{x) of excitability interval of a cell x when a number of excited 
neighbours of the cell exceeds a number of refractory neighbours. If a cell has less excited neighbours than refractory 
ones then lower boundary of the excitability interval decreases. Function £^(1, —1,0,0) increases lower boundary Oi 
and decreases upper boundary O2 of excitability interval of a cell when excited neighbours of the cells are in majority. 
If refractory cells dominate in a cell's neighbourhood then Oi decreases and O2 increases (Fig. |8|. Both functions 
decrease excitability of a cell when the cell's neighbourhood is 'over-excited' (prevalence of excited neighbours) and 
increase the cell's excitability when the cell's neighbourhood 'over-refractory' (prevalence of refractory neighbours). 
Thus they stabilise the excitation dynamics. 

Function £^( — 1,0,0,0) decreases Oi when a cell has majority of excited neighbours and increases Oi if refractory 
neighbours dominate. Function £^(—1,1,0,0) expands excitation interval (^1 decreases and O2 increases) when there 
are more excited than refractory neighbours; and contracts the excitation interval {Oi increases and O2 decreases) when 
refractory neighbours prevail (Fig. [8|. The functions increase excitability of cells with 'over-excited' neighbourhoods 
and decrease excitability of cells with 'over-refractory' neighbourhoods. Thus they destabilise the excitation dynamics. 



Finding 2 For low probabilities of initial perturbations functions stabilising excitation dynamics generate conductive 
configurations while for high probabilities of initial perturbations functions de- stabilising excitation dynamics generate 
conductive configurations. 
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FIG. 5: Fragments of configurations generated by (a) function ^(0, —1, 1, 1) from Af = (2,5) , (b) function ^(1, 1,0, 1) from 
J\f = (2, 6), (c) function E{0, 1, 1, 0) from J\f - (2, 8). 



Let £^(61,62,63,64) be a function generating conductive configurations for low probability of excitation in initial 
configurations and E\e[, 62, 63, 64) the function generating conductive configurations for high probability of excitation 
in initial configuration. Then 

£(ei, 62, 63, 64) = E\-l ' 61, -1 • 62, -1 • 63, -1 • 64) 



IV. GROWING CONDUCTIVE PATHWAYS WITH FUNCTIONS ^(1 - 100) AND ^(1000) 

Let us consider how to route conductive pathways using minimal resources. The pathways can be initiated by 
excitating the medium. A minimal seed is a pair of cells in excited state + + or ^ . The seed generates a growing 
pattern of excitation (Fig.|9|. The excitation pattern starts as a target wave (Fig.[9fc-f). At sixths step of development 
excitation ruptures inside the wave- front on its eastmost and westmost sites (Fig.|9g), and propagates along the wave- 
front boundary inside the target wave (Fig. |9]i-j). At tenth step of the seed's development the excitation reaches the 
original seed's position (Fig. lok). Thus the propagating pattern becomes 'filled' with persistent excitation activity. 



The associated development of configurations of cells in conductive state is shown in Fig. 10 The seed initiates two 
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FIG. 6: (Color online) Examples of configurations with partially connectivity clusters coloured. The clusters are bucket flooded 
with red, blue, green, magenta and yellow colours (seen as shades of grey in black and white reproduction). 



types of growing pathways: single chains and double chains of conductive states. For example, a configuration evoked 
by seed + + (Fig. [9|i) exhibits two single chains, growing west and east, and two double chains, growing north and 
south Fig. [lOj The chains growing north and south increase their lengths by one one cell each, speed 1. The chains 
growing west and east increase by one cell every other step of development, speed 1/2. 



The mechanics of formation of conductive pathways is illustrated in Fig. 11 Excited cells having configurations 



Adamatzky A., Patterns of conductivity in excitable automata with updatable intervals of excitations. Phys. Rev. E 86, 056105 (2012) [16 pages] 



10 




(c)£;(-i, 0,0,0) (d)£;(-i, 1,0,0) 

FIG. 7: (Color online) Examples of conductivity configurations generated by functions (a) ^(1,0,0,0), (b) ^(1,-1,0,0), 
(c) ^( — 1,0,0,0) and (d) ^( — 1, 1,0,0). Non-conductive cells are white. Instances of connected clusters of conductive cells are 
bucket flooded with red (c. 76 grey). 

oo— — oo oo— — oo 

o + o and o + o (north and south proximities of the excitation pattern), and configurations o + o and o + o (west 

o— oo— o oo— — oo 

and east proximities of the excitation patterns) decrease their values of Oi{x) from 2 to 1, thus increasing the cells' 
excitability. 

Finding 3 By positioning an additional cell in excited or refractory state nearby the original seed of two excited cells 
it is possible to generate a wide range of single thread wires growing into pre- determined directions. 



Adamatzky A., Patterns of conductivity in excitable automata with updatable intervals of excitations. Phys. Rev. E 86, 056105 (2012) [16 pages] 



11 



P 



Ti 



T2 



(7+ < G- 



10 



-3 



01 



02 







01 



02 



-1 



01 



02 



01 



02 



10 



-1 



-1 







01 



02 



01 



02 



-1 



01 



02 



01 



02 



FIG. 8: Mechanics of excitation interval update in functions 0,0,0), £^(1,-1,0,0), £(-1,0,0,0) and £(-1,1,0,0). cr+ is 
a number of excited neiglibours in a cell's neighbourliood and ct- is a number of refractory neighbours. 



Let us consider several examples. 

If we excite a cell above east site of the seed + + the modified seed \ ^ will develop into a growing excitation pattern 
which generates two chains of cells in conductive states. One chain grows east and another chain grows south. Space- 
time configurations of excited and refractory cells and conductive cells are shown in Fig. [12| A seed of three excited 
cells (Fig. 12 1) develops into an excitation wave fragment propagating south-east (Fig. |12^gi). The wave- fragment 
expands east and south (Fig.|12[cmo) and resembles an oxidation wave- fragment in sub-excitable Belousov-Zhabotinsky 
medium 0. Conductive wires are produced at the sites of excited cells, which have neighbourhood configurations 

o o o o — o 

° + - (north corner of the expanding wave fragment) and o + o 



(south corner) (Fig. 



12jmo). 



To produce three wires growing north, south and east we excite a neighbouring ceTTnorth-west of the seed 



(Fi 



lo pro( 
ig.|l3|- 



+ + 



The configuration 



o + 
+ + o 



of initial excitation generates a wave- fragment of excited and refractory states 



whicii travels west and expands north and south (Fig. 13 kmo). A wire growing west is produced by decreasing d\ of 

+ o - 

excited cell with neighbourhood - + o because it has more refractory neighbours (three cells) than excited neighbours 

+ o - 

o o — 

(two cells) . Wires growing north and south are produced by excited cells having neighbourhood configurations o + o 

o — o 

o — o 

and o + o . 

o o — 

By making north neighbour of eastern excited cell of the seed + + refractory we produce a single chain of conductive 
cells growing south (Fig. 14). The pattern ^ ^ (Fig. 14 i) is transformed into a particle of three excited and three 
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FIG. 9: Snapshots of a growing pattern of excitation developed from a seed of two neighbouring cells in excited state. Excited 
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FIG. 10: Snapshot of configuration of conductive cells developed from seed ++. 
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FIG. 11: Configuration of excited (•) and refractory cells (o) with highlighted cells responsible for formation of conductive 
pathways. Cells where double wire is formed are marked with arrows and single wire with arrows with disc tops. This 
configuration develops at sevenths step of the seed's evolution, Fig. (Oli. 



refractory states: 



o o — 

o o + 

— — o 

+ + o 



(Fig. 



14i)ikmoqs). The particle travels south. Its tail leaves a trace of conductive cells 



14i)ikmoqs)jlnprt) due to excited cell having two refractory neighbours: 



o + o 



and therefore decreasing its 6i 



(Fig. 
to 1. 

We have undertook an exhaustive search on head on collisions between wires growing north and south (Fig. 15). A 
wire growing south was generated by seed ^ ^ and a wire growing north by seed 1 ^ . At the beginning of each run 
seeds were placed at distance C = (/i, v) from each other, where is a number of cells along rows and v along columns 
of cellular array (Fig. p!5^). Types of collision outcomes discovered are also illustrated by schemes in in Fig. 16 ^-r. 
The following basic types are found. 

For distance C = (80, 0) the wires collide, reflect and then retract back to their origin sites. After retraction the 
continuing growing directions opposite to their original velocity vectors. Wire originally growing north grows south, 
and wire growing south grows north (Fig. [TsJ)). In the same time an additional growth seed is formed at the site of 
the collision between wires, it gives rise to two more wires growing north and south (Figs. [TsJd and [l6|^). 

In situation of one-cell horizontal shift between seeds' positions, distance C = (80, 1), both wires reflect but only 
the north growing wire (reflected south) continue growing beyond position of its original seed (Figs.jlSj^ and [l6|). If 
there is a two cell horizontal space between the seeds (Figs. [l5]i) then both wires reflect and continue their growth 
into reflected directions (Figs. 



For distance C = (80, 3) bot 



15 i and 16 n). 



1 wires reflect as in previous situations, however only one wire continues growing north 



(after passing initial position of a seed) but second wire makes 90^ turn when entering initial position of the its seed 

In situations C = (80,4) both wires just reflect without forming any 

15: and [l6b). In situations C = (80,5) 



and then grows east (Figs. 15 3 and 16 1). 



conductive bridges or patterns growing from the site of their collision (Figs, 
and C = (81,5) wires pass near each other without interacting (Figs. flStm). 

When distance between seeds is C = (81, 0) the growing wires reflect and detract back to their points of origination. 
However when they reach positions of their seeds a new pattern is formed there (Figs. 15 1). It exhibits multi-thread 
wires growing north-west, south and east from the position of northern seed, and south-west, north and east from the 
position of southern seed (Fig. 16 d). 



An E-shaped configuration of three wires growing east is formed when C = (81, 1); the seeds' site are also connected 
by a conductive wire (Figs. flSl and ^Oh). 



Finding 4 By positioning seeds 1 ^ 
wire between the seeds' locations. 



and 



+ + 



with shift C = (odd, 2) it is possible to make a stationary conductive 



For example, in situation C = (81,2) (Figs. [Ts]) the growing wires collide, make a conductive bridge and stop 
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(a) t = (b) t = (c) t = 1 (d) t = 1 




(e) t = 2 (f) t = 2 (g) t = 3 (h) i = 3 



....... OOo 



(i) i = 4 (j) t = 4 (k) i = 5 (1) t = 5 



(m) t = 6 (n) t = 6 (o) t ^ 7 (p) t = 7 



FIG. 12: Production of two wires growing east and south. For each time step t we show excitation (left pattern) and coductivity 
(right pattern) configurations. In the excitation configurations excited cells are solid black discs, refractory cells are circles. In 
the conductivity configurations conductive cells are black discs and non-conductive are grey dots. 
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FIG. 13: Three wires growing north, south and east. For each time step t we show excitation (left pattern) and coductivity 
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FIG. 14: Single wire grows south. 



(r) t = 8 (s) t = 9 (t) t = 9 

For each time step t we show excitation (left pattern) and coductivity (right pattern) 
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(a)C =(80,0) 



(b)C =(80,0) 



(c)C =(80,2) 



(d)C =(80,3) 



{e)C =(80,4) 




(f)C =(80,5) 



(g)C=(81,0) (h)C=(81,l) 
I 
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(1)C =(81,5) 



FIG. 15: Head on collision between wires growing north and south for various distances C = {h,v) between their seeds, /i is a 
number of cells between west most excited cells of seeds, and ^; is a number of cells between northmost excited cells of seeds. 
Snapshots of configurations of conductive states are taken at 270th step of iterations. The outcomes of collisions are the same 
for automata governed by functions ^(1, 0, 0, 0) and ^(1, —1, 0, 0). 



their propagation. Thus a stationary wire is formed connecting sites of southern and northern seeds (Fig. 16 :). The 
seeds' sites are connected by a wire in C = (81,3) (Figs. [Ts)^) and C = (81,4) (Figs. [Tsj) however the wire continue 
expanding north and south. 

To study outcomes of side collisions between wires we position two seeds ^ ^ and ~ ^ at distance C = (/i, v) 
(Fig. 17 1), where is a number of cells between northmost excited cells of seeds, and is a number of cells between 
northmost excited cells of seeds. The seed ^ ^ leads to formation of a single thread wire growing north. The seed 
o ^ generates a wire east. Configurations of cells in conductive states developed on 300th step after excitation of 
automaton with the seeds are shown in Fig. [17 

Colliding wires stop short of touching each other and do not propagate anymore for C = (60,40) and C = (61,43) 
(FigsJTTpk and 16 i). Growth of wire propagating north is cancelled by wire propagating east in C = (60,42) 
(Fig. [16] 



, the wires do contact each other (Figs. 17 



In condition C = (60, 43) (Figs.[T7p) two new growing wires are formed in the result of collision of a wire travelling 
north to a wire travelling east. One new wire propagates west and another wire propagates east (Fig. 



16 



Three new growing wires are formed in collision of north and east propagating wires when distance between their 
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FIG. 16: Schemes of interactions between wires. Initial positions of seeds are shown by black discs. Non-growing wires are 
shown by line segments and growing wires by arrows, (a-j) side collision, (k-r) head on collision. 



seeds is C = (60,44) (Figs. 173 andMi) and C = (61,40) (Figs. 17 1 and . One new wire propagates north and 



two new wires propagate east. The south wire propagating east is formed when originally north propagating wire 
retracts in the result of collision with originally east propagating wire and is reflected by configurations of excitation 
interval boundaries its imposed by itself. 

Almost elastic like reflection of a wire is observed in case C 



(60,45) (Figs. 17" and|l6p). A wire growing north 
collides and got cancelled by a wire growing east. In the result of impact the east growing wire is reflected and starts 



growing north-east (Fig. 16 

In situation C = (61,42) both wires retract in the result of collision. However, when they reach sites of their 
origination (where cells have already updated boundaries of excitation interval) they are transformed into extended 



17 



and 



16.) 



patterns which give rise to two new wires, both growing south (Figs. 

In situation C = (61,44) and C = (61,45) a wire growing north is retracted back to its seed's position in the result 
of collision with a wire growing east (Figs. [iTji and[l6]i). In the same time the wire growing east is reflected and 
turns north. 

Finding 5 It is possible to implement universal routing of conductive wires by positioning seeds of growing wires, the 
following operations with wires are implement able: 

• Formation of stationary wires (Fig. 

• Stopping of both growing wires (Fig.YWpj) 



Stopping of one wire by another wire without formation of conductive bridge (Fig. 16b) 
Formation of conductive circuit with one growing wire (Fig. [7^ ^ 
Formation of conductive circuit with two growing wires (Fig, fldipfikm) 
Formation of conductive circuit with three growing wires (Fig. [H 



Reflection of wires without conductive bridging (Fig. 16hno) 



Stopping of one wire and reflection of another (Fig. \Td^) 



• Co- orientation of both growing wires without formation of conductive bridge (Fig. 16g) 

• Symmetric reflection and multiplication without formation conductive bridging (Fig.^Wp). 
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(a)C =(40,0) 



(b)C =(60,40) 



(c)C =(60,42) 



(d)C =(60,43) 



^ — 



(e)C =(60,44) 



(f)C =(60,45) 



(g)C =(61,40) 



(h)C =(61,41) 



(i)C =(61,42) 



(j)C =(61,43) 



(k)C =(61,44) 



(1)C =(61,45) 



FIG. 17: Side collision between wires growing north and east for various distances C — (/i, v) between their seeds, /i is a 
number of cells between northmost excited cells of seeds, and is a number of cells between northmost excited cells of seeds. 
Snapshots of configurations of conductive states are taken at 300th step of iterations. The outcomes of collisions are the same 
for automata governed by functions E{1^ 0, 0, 0) and E{1^ —1, 0, 0). 



V. SUMMARY 



We introduced a two-dimensional excitable cellular automaton where resting cells excite depending on whether 
numbers of their excited neighbours belong to excitation intervals and boundaries of the excitation intervals are 
updated depending on ratio of excited and refractory cells in each cell's neighbourhood. We defined conductivity of 
a cell via size of its excitation interval and selected the excitation interval update functions that lead to formation of 
connected configurations of conductive cells. 

We demonstrated that by positioning elementary seeds of excitation we grow conductive wires (chains of cells in 
conductive states) and implement routing of the wires via collisions between the wires. Results presented might 
shade a light onto development of information pathways in excitable spatially extended media and contribute towards 
manufacturing of self- growing and self-organising circuits in ensembles of organic memristive polymers. 

Principle findings of the paper are following: 
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• We demonstrated that it is possible to fine tune conductivity of an excitable medium by controlling local 
dynamics of excitation. 

• Functions which stabilise excitation dynamics (where size of excitation interval increase with decrease of ex- 
citation and decreases when excitation dominates) generate fully conductive when a small number of initially 
resting cells are stimulated. 

• A point- wise initial excitation can play a seed of a growing wire, a chain of cells in conductive states; directions 
of the wire grows in pre-programmed in the configuration of initial excitation 

• The growing wires can be routed in almost arbitrary manner, dependent on positions of their seeds 

• Several wires can interact with each other by changing directions of their growth, merging in a single wire and 
co-aligning. 

We shown how to design and grow potential information pathways however we did not study how the information 
can be processed in the conductive configurations and circuits. In many cases of extended patterns are formed at the 
sites of collision between growing wires. Chances are high that these patterns can implement can implement a range 
of sensible transformations input excitation to output excitation, which could be interpreted in terms of computation. 
Computational abilities of the conductive circuits grown in excitable cellular automata will be a major topic of further 
studies. 
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